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Abstract. We study the generation of a gravitational wave (GW) background 
produced from a population of core-collapse supernovae, which form black holes 
in scenarios of structure formation of the Universe. We obtain, for example, 
that a pre-galactic population of black holes, formed at redshifts z ~ 30 — 10, 
could generate a stochastic GW background with a maximum amplitude of 
/ibg — 10 -24 in the frequency band v a ^ ~ 30 — 470Hz (considering a maximum 
efficiency of generation of GWs, namely, £gw = 7 x 10~ 4 ). In particular, we 
discuss what astrophysical information could be obtained from a positive, or even 
a negative, detection of such a GW background produced in scenarios such as 
those studied here. One of them is the possibility of obtaining the initial and final 
redshifts of the emission period from the observed spectrum of GWs. 



PACS numbers: 0430D, 9760L 

1. Introduction 

Because of the fact that gravitational waves (GWs) are produced by a large variety 
of astrophysical sources and cosmological phenomena, it is quite probable that the 
Universe is pervaded by a background of such waves. A variety of binary stars 
(ordinary, compact or combinations of them), Population III stars, phase transitions 
in the early Universe and cosmic strings are examples of sources that could produce 
such a putative background of GWs (Thorne 1987, Blair and Ju 1996, Owen et al 1998, 
Ferrari et al 1999a, 1999b, Schutz 1999, Giovannini 2000, Maggiore 2000, Schneider 
et al 2000, among others). 

In the present study we have considered the background of GWs produced from 
a Population III black hole formation. The basic arguments in favour of the existence 
of these pre-galactic black holes are the following: (a) from the Gunn-Peterson 
effect (Gunn and Peterson 1965), it is widely accepted that the Universe underwent 
a reheating (or reionization) phase between the standard recombination epoch (at 
z ~ 1000) and z > 5 as a result of the formation of the first structures of the Universe 
(see, Haiman and Loeb 1997, Loeb and Barkana 2001 for a review). However, at what 
redshift the reionization occured is still an open question (Loeb and Barkana 2001), 
although recent studies conclude that it occurred at redshifts in the range 6 < z < 30 
(Venkatesan 2000); (b) the metallicity of ~ 10~ 2 Zq found in high — z Lya forest clouds 
(Songaila and Cowie 1996, Ellison et al 2000) is consistent with a stellar population 
formed at z > 5 (Venkatesan 2000). 
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In the present paper we have adopted a stellar generation with a Salpeter initial 
mass function (IMF) as well as different stellar formation epochs. We then discuss 
what conclusions would be drawn from whether (or not) the stochastic background 
studied here is detected by forthcoming GW observatories such as LIGO and VIRGO. 

The paper is organized as follows. In section 2 we describe how to calculate 
the background of GWs produced during the formation of the stellar black holes in 
this scenario (the reader finds a more detailed description in de Araujo et al 2002), 
in section 3 we present some numerical results and the discussions, in section 4 we 
consider the detectability of this putative GW background and finally in section 5 we 
present our conclusions. 



2. The gravitational wave production from pre-galactic stars 

Before going into detail of the calculation of the background of GWs produced by 
pre-galactic stars, it is import to consider in what ways the present study differs from 
previous ones. 

Ferrari et al (1999a), for example, consider the background of GWs produced by 
the formation of black holes of galactic origin, i.e., those formed at redshifts z < 5. 
Schneider et al (2000) study the collapse of Very Massive Objects (VMO), and the 
eventual formation of very massive black holes and the GWs they generate. Fryer et 
al (2002) study, among other issues, the collapse of very massive Population III stars 
(> 300M Q ) and the GWs generated by them. 

Here, we consider a different approach to the formation of pre-galactic objects. 
We take into account the formation of Population III stellar black holes, considering 
that the progenitor stars follow a Salpeter's IMF. In this case, the progenitor stars 
have masses in the range 25 — 125Af Q . We also consider that the formation of these 
black holes occurs at high redshifts, namely, z > 10. It is worth noting that Sapeter's 
Population III stars could account for the metallicity found in high-z Lya forest clouds, 
and, at least in part, since the QSOs could also be important, for the reionization of 
the Universe. 

Let us now focus on how to calculate the background of GWs produced by the 
Population III stellar black holes we propose exist. 

The spectral energy density, the flux of GWs, received on Earth, F v , in 
erg cm -2 s _1 Hz -1 , reads (see, e.g. Douglass and Braginsky 1979, Hils et al 1990) 

p _ c3s h^obs m 

lbnG 

where w bs = 27W bs, with v Q \> s being the GW frequency (Hz) observed on Earth, c is 
the velocity of light, G is the gravitational constant and ^/sh is the strain amplitude 
of the GW (Hz" 1 / 2 ). 

The stochastic GW background produced by gravitational collapses that lead to 
black holes would have a spectral density of the flux of GWs and strain amplitude 
also related to equation (|l|) . The strain amplitude at a given frequency, at the present 
time, is the contribution of black holes with different masses at different redshifts. 
Thus, the ensemble of black holes formed produces a background whose characteristic 
strain amplitude at the present time is \^s h . 
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On the other hand, the spectral density of the flux can be written as (Ferrari et 
al 1999a) 

F v = / / fu(vohs)dRBH(m, z) (2) 

J Z c { J m m i n 

where f u (v bs) is the energy flux per unit of frequency (in erg cm -2 Hz -1 ) produced 
due to the formation of a unique black hole and cLRbh is the differential rate of black 
hole formation. 

The above equation takes into account the contribution of different masses that 
collapse to form black holes occurring between redshifts z c \ and z c { (beginning and end 
of the star formation phase, respectively) that produce a signal at the same frequency 
fobs- On the other hand, we can write /^(fobs) (Carr 1980) as 

/*Oob B ) = ^Q h BH ( 3 ) 

where /ibh is the dimensionless amplitude produced by the collapse, to a black hole, 
of a given star with mass m that generates at the present time a signal with frequency 
z/ D bs- Then, the resulting equation for the spectral density of the flux is 

7TC 3 f 

FV = 2G J ^H^BH. (4) 

From the above equations, we obtain for the strain amplitude 

Sh = / h% n dRBH- (5) 

^obs J 

Then, the dimensionless amplitude of the background of GWs, /isd reads 

h% G = — I ft.| H di?BH (6) 

^obs J 

(see, de Araujo et al 2000 for details) where f bsj ^bh arid oIRbk are defined below. 

The dimensionless amplitude ft-BH produced by the collapse of a star, or star 
cluster, to form a black hole is (Thorne 1987) 

where £qw is the efficiency of generation of GWs, M r is the remnant black hole mass 
and (iL is the luminosity distance to the source. 

It is worth mentioning that equation (0) refers to the black hole 'ringing', which 
has to do with the de-excitation of the black hole quasi-normal modes. Note also that 
£gw ^ fl4 (see, e.g. Stark and Piran 1986), where 'a' is the the dimensionless angular 
momentum. Thus, greater the GW efficiency, greater the dimensionless angular 
momentum. 

We assume that the progenitor masses of the black holes range from M = 
25 — 125M Q (see Timmes et al 1995, Woosley and Timmes 1996). The remnant 
and the progenitor masses are related to M r = aM and, we assume a — 0.1 (see, e.g. 
Ferrari et al 1999a). 

The collapse of a star to a black hole produces a signal with an observed frequency 
v ohs at the Earth (Thorne 1987) 

^ obs ~ 1.3 x H^Hzf^Vl + z)- 1 (8) 
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where the factor (1 + z)~ x takes into account the redshift effect on the emission 
frequency. 

The differential rate of black hole formation di?BH reads 

oIRbh — p*(z)-^-(j)(m)dmdz (9) 

where dV is the co-moving volume element, <f>{m) is the stellar initial mass function 
(IMF) and is the star formation rate (SFR) density. 

The SFR density can be related to the reionization of the Universe. The 
amount of baryons necessary to participate in early star formation, to account for 
the reionization, would amount to a small fraction, of all baryons of the Universe 
(see, e.g. Loeb and Barkana 2001). We then assume that 

= % = < 10 > 

where the term in brackets represents the stellar mass density at redshift z, with p c 
the present critical density, and fi* the stellar density parameter. The latter can be 
written as a fraction of the baryonic density parameter, namely, f2+ = /*Ob, which 
we assume to be independent of the redshift. 

From the above equations, we obtain for the dimensionless amplitude 



2 (7.4 x 10- 2 °q) 2 £GW 

t'obs 

where m m i n = 25Mq, m u = 125Mq and z c \ (z c f) is the beginning (end) of the black 
hole formation phase. Equation ([ll]) is computed for each observed frequency. Also, 
looking at equation (pT|), one notes that to integrate it, one needs to choose the IMF, 
the cosmological parameters and set values for the following parameters: z C i, z c [, a, 
£gw, /+■ In the next section we present the numerical results and discussions. 



3. Numerical results and discussions 



To evaluate the background of GWs produced by the formation of the Population 
III black holes, it is necessary to know in which redshifts they began and finished 
formation. This is a very difficult question to answer, since it involves knowledge of 
the role of the negative and positive feedbacks of star formation which are regulated 
by cooling and injection of energy processes. 

Should the stochastic background of GWs studied here be significantly produced 
and detected at a reasonable confidence level, the present study can be used to obtain 
the redshift range where the Population III black holes were formed. We refer the 
reader to the paper by de Araujo et al (2002) for further discussions. In figure 1 an 
example is given of how one could get z c \ and z c f from the curve /ibg versus v bs- 
Knowing the frequency band v min — ^max detected from a cosmological source, and 
using equation (^), one can obtain both z ci and z c { from figure 1. Thus, these redshifts 
are therefore observable. Note that we have assumed as did Ferrari et al (1999a) that 
a is a constant (a = 0.1). 
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Figure 1. Example of how one could obtain the initial (final) redshift z c \ (« c f) 
of the GW emission period from /ibg versus f bs- We have adopted a = 0.1. 



Stars start forming at different redshifts, creating ionized bubbles (Stromgren 
spheres) around themselves, which expand into the intergalactic medium (IGM), at a 
rate dictated by the source luminosity and the background IGM density (Loeb and 
Barkana 2001). The reionization is complete when the bubbles overlap to fill the entire 
Universe. Thus the epoch of reionization is not the epoch of star formation. There is 
a non-negligible time span between them. Here, we have chosen different formation 
epochs to see their influence on the putative background of GWs and also to see if it 
could be detected by the forthcoming GW antennas. 

To calculate /ibg we adopted the standard Salpeter IMF. For £gw, the efficiency 
of production of GWs, whose distribution function is unknown, we have parametrized 
our results in terms of its maximum value, namely, £GW max — 7 x 10~ 4 . This figure 
is obtained from studies by Stark and Piran (1986) who simulated the axisymmetric 
collapse of a rotating star to a black hole. 

To calculate /ibg we still need to know fi*, which has a key role in the definition of 
the SFR density. From different studies one can conclude that a few per cent, may be 
up to ~ 10%, of the baryons must be condensed into stars in order for the reionization 
of the Universe to take place (see, e.g. Venkatesan 2000). Here we have set the value 
of f2* in such a way that it amounts to 1% of all baryons (our fiducial value). 

Looking at equation (jll]), one could think it would depend critically on the 
cosmological parameters Hq, f2e, ^dm (the density parameter for the dark matter) 
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Table 1. The redshifts of collapse for our models and the corresponding GW 
frequency bands. The cosmological parameter h^ 00 QB = 0.019 (see the text), 
a = 0.1, /* = 0.01 (our fiducial value) and the standard IMF arc adopted. 



Model 




z c f 


Ai/ obs (Hz) 


A 


20 


10 


50-470 


B 


30 


20 


34-250 


C 


40 


30 


25-170 


D 


50 


40 


20-130 


E 


30 


10 


34-470 


F 


40 


10 


25-470 


G 


50 


10 


20-470 



and ft a (the density parameter associated with the cosmological constant). But our 
results show that /ibg depends only on Hq and Qb, the Hubble parameter and the 
baryonic density parameter, respectively. 

The quantity /i^qq^b = 0.019 ±0.0024 (where h\oo is the Hubble parameter given 
in terms of 100 km s _1 Mpc -1 ) is adopted in the present model. This figure is obtained 
from Big Bang nucleosynthesis studies (see, e.g. Buries et al 1999). 

In table 1 we present the redshift band, z c ; and z c {, for the models studied and 
the corresponding GW frequency bands. For the cosmological parameters, we have 
adopted hioo — 0.65, CIm = 0.3, f2e = 0.045 and Oa = 0.7. We have also adopted 
a = 0.1, /* = 0.01 and the standard IMF. 

It is worth noting that no structure formation model has been used to find the 
black holes formation epoch; instead, we have simply chosen the values of z to see 
whether it is possible to obtain detectable GW signals. In the next section, it will 
be seen that unless £gw is negligible, the background of GWs we propose here can 
be detected. Our choices, however, can be understood as follows. The greater the 
redshift formation, the more power the masses related to the Population III objects 
have. Thus from our models A to D, our model D (A) has more (less) power when 
compared to the others. The models E, F and G would mean a more extended star 
formation epoch, which means that the feedback processes of star formation are such 
that they allow a more extended star formation epoch when compared to the models 
B, C and D, respectively. 

Note that the reionization epoch occurred at lower redshifts as compared to the 
first stars formation redshifts. Loeb and Barkana (2001) found, for example, that if 
the stars were formed at z ~ 10 — 30, with standard IMF, they could have reionized 
the Universe at redshift z ~ 6. Our models A, B and E, for example, could account 
for such a reionization redshift. 

If the process of structure formation of the Universe and the consequent star 
formation were well known, one could obtain the redshift formation epoch of the first 
stars. On the other hand, if the background of GWs really exists and is detected, one 
can obtain information about the formation epoch of the first stars. 

A relevant question is whether the background we study here is continuous or 
not. The duty cycle indicates if the collective effect of the bursts of GWs generated 
during the collapse of a progenitor star generates a continuous background. For all 
the models studied here the duty cycle is ^ 1 (see de Araujo et al 2002 for details). 

We find, for example, that the formation of a Population (III) of black holes, 
in the model D, could generate a stochastic background of GWs with amplitude 
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Table 2. For the models in table 1, we present the SNR for pairs of LIGO I, II 
and III ('first', 'enhanced' and 'advanced', respectively) observatories for 1 year of 
observation. Note that an efficiency of generation £GW max = 7 X 10~ 4 is assumed. 



Model 




SNR 




LIGO I 


LIGO II 


LIGO III 


A 


8.3 x 1CT 3 


1.6 


6.6 


B 


8.5 x 1CT 3 


2.3 


26 


C 


8.7 x 1CT 3 


2.7 


47 


D 


8.1 x icr 3 


2.5 


51 


E 


2.7 x icr 3 


5.7 


37 


F 


5.0 x icr 3 


12 


120 


G 


7.7 x icr 2 


21 


260 



/ibg — (0.8— 2)xl0~ 24 and a corresponding closure density of SIqw — (0.7— 1.4) x 10 -8 , 
at the frequency band j/ b s — 20 — 130 Hz (assuming an efficiency of generation 
£qw — 7 x 10~ 4 , the maximum one). 

In other paper to appear elsewhere (de Araujo et al 2002), we study in detail how 
the variations of the several parameters modify our results. 

4. Detectability of the background of gravitational waves 

The background predicted in the present study cannot be detected by single 
forthcoming interferometric detectors, such as VIRGO and LIGO (even by advanced 
ones). However, it is possible to correlate the signal of two or more detectors to detect 
the background that we propose to exist. 

To assess the detectability of a GW signal, one must evaluate the signal-to-noise 
ratio (SNR), which for a pair of interferometers is given by (see, e.g. Flanagan 1993) 



SNR 2 = 



(S)*f* 



(12) 



where is the spectral noise density, T is the integration time and ^{v) is the 
overlap reduction function, which depends on the relative positions and orientations 
of the two interferometers. The closure energy density is given by 



1 dp GW 47T 2 2 

/9 C dlog^obs 3Hf 



Here we consider, in particular, the LIGO interferometers, and their spectral noise 
densities have been taken from a paper by Owen et al (1998). 

In table 2 we present the SNR for 1 year of observation with a = 0.1, SIb^-ioo = 
0.019, /* = 0.01 and £ G w max = 7 x 10~ 4 for the models in table 1, for the three LIGO 
interferometer configurations. 

Note that for the 'initial' LIGO (LIGO I), there is no hope of detecting the 
background of GWs we propose here. For the 'enhanced' LIGO (LIGO II) there is 
some possibility of detecting the background, since SNR > 1, if e GW is around the 
maximum value. Even if the LIGO II interferometers cannot detect such a background, 
it will be possible to constrain the efficiency of GW production. 
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The prospect of detection with the 'advanced' LIGO (LIGO III) interferometers 
is much more optimistic, since the SNR for almost all models is significantly greater 
than unity. Only if the value of £qw were significantly lower than the maximum 
value would the detection not be possible. In fact, the signal-to-noise ratio is critically 
dependent on this parameter. 

Note that the larger the star formation rcdshift band, the greater the SNR. 
Secondly, the earlier the star formation, the greater the SNR. It is worth recalling 
that, if one can obtain the curve '/ibg versus fobs' and the value of a is known, one 
can find the redshift of star formation. 

5. Conclusions 

We have shown that a background of GWs is produced from Population III black 
hole formation at high rcdshift. This background can in principle be detected by a 
pair of LIGO II (or more probably by a pair of LIGO III) interferometers. However, a 
relevant question should be considered: what astrophysical information can one obtain 
whether or not such a putative background is detected? 

First, let us consider a non-detection of the GW background. The critical 
parameter to be constrained here is ecw- A non-detection would mean that the 
efficiency of GWs during the formation of black holes is not high enough. Another 
possibility is that the first generation of stars is such that the black holes formed had 
masses > 100M Q , and should they form at z > 10 the GW frequency band would be 
out of the LIGO frequency band. 

Secondly, a detection of the background with a significant SNR would permit 
us to obtain the curve of /ibg versus v b s - From it, one can constrain a and the 
rcdshift formation epoch; and for a given IMF and fie ft 2 , one can also constrain the 
values of /+ and £gw- On the other hand, using the curve of ft.BG versus ^ bs and in 
addition other astrophysical data, say CBR data, models of structure formation and 
reionization of the Universe, constraint on £gw can also be imposed. We refer the 
reader to the paper by de Araujo et al (2002) for further discussions. 

It is worth mentioning that a significant amount of GWs can also be produced 
during the formation of neutron stars and if such stars are r-mode unstable (Andcrsson 
1998). We leave these issues for other studies to appear elsewhere. 
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